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ABSTRACT 

The use of the tensor virial theorem (TVT) as a diagnostic of anisotropic velocity 
distributions in galaxies is revisited. The TVT provides a rigorous global link between 
velocity anisotropy, rotation and shape, but the quantities appearing in it are not easily 
estimated observationally. Traditionally use has been made of a centrally averaged 
velocity dispersion and the peak rotation velocity. Although this procedure cannot 
be rigorously justified, tests on model galaxies show that it works surprisingly well. 
With the advent of integral-field spectroscopy it is now possible to establish a rigorous 
connection between the TVT and observations. The TVT is reformulated in terms of 
sky-averages, and the new formulation is tested on model galaxies. 
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1 INTRODUCTION 

Thirty years ago our understanding of elliptical galaxies was 
revolutionized by the discovery that most giant elliptical 
galax i es are not flattened by rotat i on l|Bertola fc Capaccioil 
Il977t llllingworthl Il977t iBinnevI Il978h . Subsequently it 
emerged that a galaxy's peak rotation speed is correlated 
with its luminosity, and the cuspiness and diskiness of its lu- 
minosity density: less luminous galaxies tend to have cuspier 
central density profiles, disky rather than boxy isophotes, 
and a higher d egree of rotat i onal flattening than more h imi- 
nous galaxies <BendeJll983 iKormendv fc Benderlll99r3) . 

It is generally recognized that these correlations must 
be clues to how elliptical galaxies formed. A promising con- 
jecture is that the current configuration of a boxy galaxy 
was largely established by a violent merger of galaxies of 
comparable mass, while cuspy galaxies were configured by a 
succession of minor mergers wi th companions sub stantially 
less massive than themselv es (iNaab. Burkert fc Hernauisd 
ll999l:lNaab fc BurkertfeOOSft . To assess the truth of this and 
any competing conjecture, we clearly need to characterize 
as fully as possible the degree of rotational flattening of any 
given galaxy, and to understand the significance of rotational 
flattening for the system's internal dynamics. 

Ideally one would explore the connections between disk- 
iness, cuspiness and anisotropy by experimenting with a 
series of semi-analytic galaxy models that wo uld be gen- 
eralizations of popular sph erical models (e.g.. lKin3ll966t 
|Jaffdll983l: iHernauisd 11990). Unfortunately, even now the 
repertoire of semi-analytic flattened models is very sparse, 
and all of these models are unrealistic in one significant 
respect or another. Moreover, thirty years ago, when the 
importance of anisotropy was first realized, even fewer 
semi-analytic models were available, and our theoretical 



understanding has relied heavily on (a) N-body models 
j |Binnevlll976t [Aarseth fc Binnevlll978l : [Bjirngs^^TIernrjuisd 
I l992ft . (b) Schwarzschild modelling JSchwarzschildl ll97St 
iLevison fc Richstonel lT980l) . and (c) the tensor virial the- 
orem feinnevlll97gFl . 

The tensor virial theorem (TVT) provides a powerful 
general framework within which to discuss the connection 
between rotation and flattening, but it is not straightfor- 
ward to connect the quantities that appear in it to observ- 
able quantities. Until recently it was in principle impossible 
to make this connection rigorously because the TVT op- 
erates at a global level, while galaxy kinematics could be 
probed only along a limited number of slits, that rarely ex- 
tended very far out in the galaxy, especially perpendicular 
to the apparent major axis. With the advent of integral- 
field spectroscop y jBac on et al.ll2002t Ide Zeeuw et aljl2002l : 
iKelz. Roth fc Beckerll2005t this problem is substantially al- 
leviated, and it is now time to reassess how the TVT is used 
to interpret kinematic data. This task is the primary goal of 
this paper. In the process we make a critical reassessment 
of how traditional long-slit data are interpreted. 



2 THE TENSOR VIRIAL THEOREM 

The TVT states that in any equilibrium ste llar system 
<Chandrasekhadll969l iBinnev fc Tremainelll987h 

2K + W = 0, (1) 

where the kinetic- and potential-energy tensors are 



Ki 



d w ViVjf(w) 
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ir, 



' 3 a$ 



(2) 



Here w = (x, v) is the vector of phase-space coordinates, 
/ is the system's distribution function, p(x) = Jd 3 v/(w) 
is the density, and ( I > (x) is the gravitational potential. K 
is customarily decomposed into contributions from ordered 
and random motion 



K 

where 

Ti 1 



T + fn, 



(3) 



§ / d 3 xw., 



n, 



d w(vi -Vi)(vj -Vj)f(-w). 



(4) 



Here a bar over any quantity denotes an average over veloc- 
ity space: 

^(x) = ^y"d 3 v^/(w). (5) 
We use this notation to define the velocity-dispersion tensor 



Vi)(Vj 



(6) 



In the case of a flattened, axisymmetric galaxy, we ex- 
pect the principal axes of the tensors above to coincide with 
the symmetry axis, which we label the z axis, and any two 
perpendicular axes - the x and y axes. Then II can be char- 
acterized by two numbers: a\ and the global anisotropy pa- 
rameter which are defined such that 



n zz = (1 



S)Ma 2 Q , 



(7) 



where M is the galaxy's mass. If we assume that the galaxy's 
only streaming motion is rotation, we may characterize T 
by a single number vq through 



\Mvt 



(8) 



where the factor j is chosen so that the total ordered kinetic 
energy is T xx + T yy — ^Mv 2 . With thi s notation it i s easy 
to show that the TVT Q implies that feinnevlll97sl) 

^|=2(l-5)^-2. (9) 



In the case that the galaxy's surfaces of constant density are 
all similar spheroids of axis ratio 1 — e, the ratio W xx /W zz 
on the right of this equation is a function of e and indepen - 
dent on the galaxy's radial density profile dRobertslll962l) . 
so this equation provides a connection between the rota- 
tion rate vo, the mean-square velocity dispersion parallel to 
the equatorial plane a 2 , and the global anisotropy param- 
eter S. If «o and 00 could be estimated from observational 
measurements of the line-of-sight mean velocity and velocity 
dispersion, equation @ would enable one to determine the 
anisotropy parameter 5 from observational data. 



3 APPLICATION TO TRADITIONAL DATA 



Following lirlingworthl Jl977l) and lDavies et al.Hl98ot . a has 
often been identified with 00.5, the mean velocity dispersion 
interior to ~ ^Rc, while the peak line-of-sight streaming ve- 
locity Umax is identified with jttvo- There is no rigorous ba- 
sis for either identification since both the velocity dispersion 
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Figure 1. Squares: the ratio of peak line-of-sight streaming ve- 
locity to line-of-sight velocity dispersion for two Evans models. 
Triangles: the same data for three Rowley models. The point at 
the upper right of each chain is for an edge-on model, and suc- 
cessive points on the chain show the effect of reducing the incli- 
nation by 10°. The dashed curve shows the relation J^J between 
u maI = 77TVo and (to for <5 = 0. The circles show observed values 
of i)max/"0.5 for spheroids less luminous than Mg = —20.5. 



and the streaming velocity are expected to depend on the 
location of the line-of-sight. The factor j-k is motivated by 
the consideration that if the luminosity density j(R, z) scales 
with galactocentric distance (R 2 + z 2 )~ 3 ^ 2 , and the stream- 
ing velocity is independent of R, then when the galaxy is 
viewed edge-on, the line-of-sight streaming velocity will be 
a factor 7r/4 smaller t han the three-dimensional streaming 
velocity <Binnevlll97sh . 

We can test the reliability of these traditional identifi- 
cations by applying them to models that have analytic dis- 
tribution functions of the form f(E,L z ) and self-consistent 
gravitational potentials, where E and L z are respectively 
the stellar energy and angular momentum about the sym- 
metry axis. Fig.Qshows data for two types of semi-analytic 
model. 



3.1 Evans models 



lEvansl il994f) presented a family of models, called 'power- 
law models', in which the potential varies as a power of 
the spheroidal variable m = R 2 + R 2 + z 2 /q\ . Here R c is 
a constant that sets the linear scale of the model, while qg> 
determines the model's flattening. The index y in the power- 
law relation $ oc (R 2 /m 2 ) v determines the asymptotic slope 
of the circular-speed curve - the models shown in Fig. 
both have y — 0.09, which corresponds to a gently falling 
asymptotic circular-speed curve. 

The part f+(E,L z ) of the distribution function that is 
even in L z is uniquely determined by Poisson's equation, 
and we define the odd function f-(E,L z ) to be equal to 
f+(E,L z ) for L z > 0. The complete distribution function is 
taken to be / = /+ + 7/- . With this choice of distribution 
function, the streaming velocity at any point in the model 
is proportional to 7 and we choose 7 as the average of the 
individual values that would make — or at a grid of 
points in the equatorial plane. That is, 7 is chosen to make 
the model as near as possible an isotropic rotator. 
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The squares in Fig. show centrally averaged line-of- 
sight velocity dispersion and v m&K for Evans models with 
gj. = 0.85 and g$ = 0.9. The value of e is for the isophote 
with semi-major axis length 15R C , and cro.s is the mean value 
of the line-of-sight velocity dispersion on the major axis out 
to 15-Rc- The filled squares are for a model seen edge-on, 
while successive open squares show data for the models seen 
at inclination i = 80° , 70° , . . . . 

3.2 Rowley models 

iRowlevI il988t) developed models of flattened stellar systems 
to use as model bulges. These systems are defined by their 
distribution function 



for x > 0, 
otherwise, 



where 



X = Xo - E + ujL z - \L 2 z /r 2 a , 



(10) 



(11) 



with E stellar energy, L z angular momentum about the sym- 
metry axis, and a, \o, w and r a being constants. It is easy 
to show that x can be rewritten as 



X = 



*-K-|(i + ; |-) («*- T 



ujR 



(12) 



(13) 



+ R 2 /rl. 

where v m is the speed in the meridional plane, and 
lo 2 R 2 

X ° 2(1 + B?/riy 

From equations 11011 and 1 1 2t we can identify the streaming 
velocity 

«* - mm* (14) 

and observe that the distribution of stellar velocities is a 
truncated bi-axial Gaussian around v^. We see also that the 
density vanishes for <]/ < 0. We define the tidal radius Rt 
as the radius at which this inequality is first satisfied in the 
equatorial plane. Rowley models are characterized by the 
dimensionless spin parameter ior a /a 2 and the dimensionless 
potential difference A = [*(0, 0) - V(R t , 0)]/<r 2 . 

Each filled triangle in Fig. shows the ratio of u nB to 
centrally-averaged velocity dispersion versus the ellipticity 
of the isophote at 3R C for a Rowley model seen edge-on, 
while the open triangles show the corresponding data for the 
models seen at inclinations i = 80°, 70°, .... All models have 
A = 4, and u>r a /cr — 0.9, 1.35, 1.69 for the three models. 

We see that at a given ellipticity, Rowley models rotate 
more rapidly than Evans models. This phenomenon reflects 
the low ratio of to gr that follows from equation 1121 . 
When a razor-thin disk is added to a model, the flattening 
of the system is increased while its rotation stays the same, 
and for realistic parameters the system can move into the 
band in the (v/a, e) plane that is occupied by bulges iRowlev| 
Il988h . 

The dashed curve in Fig. Q shows the relation be- 
tween j-KVo/ao and e that is predicted by equation ©. The 
curve runs close to the location of edge-on Evans models, 
which are very nearly isotropic rotators, while the Row- 
ley models, which are slightly anisotropic in the sense that 
< a \ = °zj lie above the dashed curve. Thus Fig. [T] 



validates the traditional use of the TVT to draw the trajec- 
tory of isotropic rotators in the (v/a,e) plane. The models 
demonstrate that tilting a galaxy away from edge-on ori- 
entation tends to push its representative point above the 
dashed curve. Hence galaxies that lie well below this curve, 
as luminous elliptical galaxies nearly all do, must certainly 
be anisotropic in the sense of having larger pressure parallel 
to the equatorial plane than perpendicular to it. 

In Fig. circles show ratios of peak rotation speed to 
velocity dispersion averaged within ^R e for real spheroidal 
systems less luminous than Mb = —20 .5 (for Ho = 
50kms _1 Mpc" 1 ) from iDavies et al.l Jl983l) . These circles 
straddle the dashed curve, with rather more points below 
than above it, suggesting that these systems are at most 
mildly anisotropic. 



4 APPLICATION TO MODERN DATA 

Observers measure the line-of-sight velocity distribution 
(LOSVD) 



F(x_ Ll ui|) = 



E(xx) 



da; 



d 2 v x /(w), 



where 

E = I das || p(x) 



(15) 



(16) 



is the surface density and the subscripts || and _L respectively 
denote components of vectors parallel and perpendicular to 
the line of sight. They quantify F by the moments 



a 2 (xx) = / dwy («|| - «n) 2 F(xj_,«||) 



(17) 



3||( X ^)= / dv \\ S( W I|) F ( X -L> U I|) 

= h / dx \l / d3v 9( v \\ )/( w ) 



and v\\ (x_l), where we have introduced the notation that for 
any function g(vn), 



(18) 



Integral-field spectr ographs such as SAURON 
ide Zeeuw et alj 12002?) enable one to map vy and a 2 
over a significant part of the galaxy image. From such 
maps it is possible to measure quantities that are more 
directly related to the TVT than the quantities fj max and 
o"o. 5 defined above. 

First we note that 



2 2 ~2 

CTll = v„ — V\\ . 



(19) 



Using this relation to integrate Su 2 over the sky, we obtain 



d 2 x ± (of +£?)£= / d 3 x / d 3 v«jj/(w) 



= 2s • K ■ s, 



(20) 



where s is the unit vector parallel to the line of sight. 

For s parallel to the as axis we can write this equation 
more compactly as 



Of) + <«!>) M = 2K **> 



(21) 
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Figure 2. Modified TVT. The full curves show the relation be- 
tween (t>|)/(o"|) established by equation 1261 when a = 0, while 
the dotted curves show the relation when a = 0.2. 



where we have introduced the notation that with q(x.±) an 
arbitrary function on the sky, its sky-average is 



<9> = M / d " X±gE ' 



(22) 



Equation 12 H for K xx enables us to derive from the xx and 
zz components of the TVT 



M 



(23) 



At any point x in the galaxy let tt(x) be the difference 
between the component V« of the streaming velocity parallel 
to the line of sight and the mean velocity for that line of 
sight, vu. With this definition, u = Vii — v\\, it is easy to 
show that 



San 



dxn (a xx + w 2 )p(x), 



(24) 



where a xx is defined by equation JBJ. Integrating this ex- 
pression over the sky, we obtain 



!!« = M(a?i) - / d 3 xu 2 (x)p(x) 



(25) 



Finally we use this equation and equation (J) to eliminate 
II 22 from equation 12311 . and obtain after some rearrange- 
ment 



M> _ 

where 



(1 - S)W XX /W ZZ - 1 
a(l-*)W™/W„ + l' 



d xa p( x ) 



(26) 



(27) 



In this form of the TVT, the left side contains only observa- 
tionally accessible quantities, while the right side contains 
the standard global anisotropy parameter 5, the shape pa- 
rameter W xx /W zz , and a new dimensionless parameter a, 
which quantifies the contribution of streaming motion to 
the line-of-sight velocity dispersion. 

The value of a depends on the shape of the galaxy's 
intrinsic rotation curve, and its radial luminosity profile, 
but not on the amplitude of the rotation curve. One may 



straightforwardly show that it vanishes in the case of solid- 
body rotation, and increases with the shear of the stellar 
flow. For a fixed rotation curve, it is larger for flatter slopes 
of the radial luminosity density profile. A numerical calcula- 
tion for th e case of a galaxy in which the luminosity density 
follows the lHernauisd (jl990T ) profile shows that a = 0.131 if 
«0 = constant. 

Fig. [5] shows v/a as a function of e from equation 1261 
for several values of 8 and two values of a: a = (full curves) 
and a = 0.2 (dotted curves). When a — the right side of 
equation 126H becomes precisely half of the right side of the 
classical result JSJ. The factor half arises because in equation 
© wo is a three-dimensional streaming velocity, whereas in 

equation 126H \J (vf\) i s a projected velocity. 

Probably equation (12(il is most usefully recast as an 
expression for the anisotropy parameter S as a function of 



and {v >/(a >: 



6 = 1- 



l + <«f>/K 



V_ fW zz \ 
tf) \ w xx ) 



(28) 



4.1 Practical considerations 

Three issues arise when using observational data to evaluate 
the right side of equation 1281 . We cannot directly measure 
a, but we can infer its value with reasonable accuracy from 
the shape of the measured rotation curve. In the tests below, 
I simply use the values directly calculated from the models. 
The second issue is that our kinematic data do not extend 
to arbitrary radii, so it is not practicable to evaluate the 
means and (fn). In practice our averages must be con- 
fined to some inner region and we must investigate the seri- 
ousness of the error incurred by this confinement. The final 
problem is how to determine the required ratio W zz /W xx of 
components of the potential-energy tens or. In principle this 
may be done from the photometry (e.g. iBinnev fc Strimpell 
1978). A much simpler alternative is to assume that the iso- 
density surfaces have a constant ellipticity e, which allows 
W zz /W xx to be determined independently of the radial den- 
sity profile. In most of the tests below I have used this crude 
approach, with e taken to be a weighted average e of the 
ellipticities of individual isophotes. After some experimen- 
tation, the weighting scheme adopted was 



/ dRR 2 T,(R,0)e(R) 
JdRR 2 X(R,0) ' 



(29) 



where E(7i, 0) denotes the surface brightness distance R 
down the apparent major axis. 

Insight into the uncertainties that will be encountered 
when using real data is provided by applying equation 128> 
to pseudo-data obtained by projecting Evans and Rowley 
models. 

The upper panel of Fig. |H] shows the projected rota- 
tion rate and ellipticity of an Evans model with y = 0.25 
and q$ = O.9. 1 Because the rotation curve rises very steeply 
near the centre (in three dimensions v does not tend to zero 
as the origin is approached), a ~ 0.2 is relatively large: the 



1 The final integral in the definition of Wij diverges for Evans 
models with y < 0.25. 
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Figure 3. Upper panel: projected rotation speed and cllipticity 
for an Evans model. Lower panel: the parameter a denned by 
equation 1271 (dotted curve) and the global anisotropy parame- 
ter <5 estimated from equation 1281 using only data from within 
the circle of radius R; at edge-on orientation (full curve) and for 
inclination i = 30° (short-dashed curve). The long-dashed curve 
shows values of 5 recovered at edge-on orientation when the true 
ratio W zz /W x:l: is used in 1281 instead of an estimate based on 
the assumption that all isodensity surfaces are similar spheroids. 



dotted curve in the lower panel of Fig. [3]shows as a function 
of radius R the value of a obtained by carrying the integrals 
involved in its definition l|27|l only out to R. The full curve 
in the lower panel shows as a function of R the values of 5 
yielded by equation 1281 when the model is seen edge-on and 
W Z z/W xx is evaluated under the assumption of constant el- 
lipticity, e. Small positive values of 8 < 0.075 are obtained, 
and outside the core (where a becomes large) <5 rises slowly 
with R because e is rising, while the streaming velocity is 
very nearly constant. The long-dashed curve shows the val- 
ues of 8 obtained when W zz /W xx is set equal to the value it 
assumes at large radii. The decline of this curve reflects the 
outward-increasing nature of e. The fact that at R ~ 50R C 
this curve has fallen to small negative 8 ~ —0.05 suggests 
that the procedure intended to make the model an isotropic 
rotator was not wholly successful. The offset between the 
end points of the full and long-dashed curves gives an indi- 
cation of the error in S that is inherent in adopting a constant 
ellipticity. 

Fig. [IJshows equivalent data for the Rowley model that 
has A = 4 and lot a /a — 1.69. Because this model has a sub- 
stantial region of near solid-body rotation, a <. 0.03 is much 
smaller than in the case of the Evans model. When the model 
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Figure 4. As Fig. [3] but for a Rowley model. 



is viewed edge-on, 8 falls by 0.35i?t to 8 = —0.075 and from 
there to the edge of the model oscillates in a small range 
around —0.097. The end-point of the long-dashed curve, 
which is based on the exact ratio W Z z/W xx , lies within this 
range because at most radii the azimuthal velocity disper- 
sion is smaller than the other two principal dispersions. 

The short-dashed curves in the lower panels of Figs [3] 
and 0] show that remarkably small values of S are obtained 
when equation 1281 is applied to data obtained from the 
models at inclination i = 30°. Since equation 128H was de- 
rived for edge-on inclination, we have no guarantee that a re- 
alistic value of 8 will be obtained at small inclinations. More- 
over, in Fig. both models move away from the isotropic- 
rotator line towards negative 8 as the inclination is reduced. 
Hence it comes as a pleasant surprise to find that equation 
128t works well at small inclinations. 



5 CONCLUSIONS 

The degree of anisotropy in the velocity distribution of an el- 
liptical galaxy is probably an important clue to the manner 
in which the galaxy formed. Traditionally anisotropy has 
been estimated from the ratio w ma x/o"o.5 of the maximum 
measured rotation velocity to a central mean of the velocity 
dispersion. By a distinctly ad-hoc scaling, the tensor virial 
theorem has been used to draw a curve on the (v/a, e) plane 
for isotropic rotators, and galaxies that lie significantly be- 
low this curve have been deemed to have anisotropic veloc- 
ity distributions. Fig. ^validates this procedure by showing 
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the models with nearly isotropic velocity distributions are 
placed near the curve of isotropic rotators. 

A much more rigorous procedure is possible now that 
one can map rotation speed and velocity dispersion over 
a substantial fraction of a galaxy's image. One replaces 
"max/cro.5 with the ratio of the sky-averaged 

squared rotation velocity to the squared velocity disper- 
sion. Equation 12611 rigorously relates this average to the 
global anisotropy parameter 8, a shape-dependent ratio of 
potential-energy tensor components, and the dimensionless 
quantity a that measures the degree of shear in the stel- 
lar streaming velocity. Alternatively, equation 12811 expresses 
the anisotropy parameter as a function of quantities that are 
either directly observable ((«y)/(o"y) and W zz /W xx ) or can 
be estimated to sufficient accuracy from the data (a). 

An obvious problem with the rigorous approach is that 
the sky-averages that appear in the TVT extend over the 
whole image, while data are available only within some limit- 
ing radius R. Figs|3]and2]suggest that 5 can be estimated to 
an uncertainty ~ 0.05 from data that extend out to > 15R C 
in a galaxy that resembles an Evans model (in which the 
ellipticity reaches a plateau around 10-R C ), or >0.357?t in 
a galaxy like a Rowley model (in which e roughly linearly 
to the edge of the system). It is unfortunate that when we 
confine ourselves to semi-analytic models, we are able to 
investigate the impact of restricted sky coverage only with 
models that have 5 ~ 0. Tests on systems with i > could 
be carried out with N-body models. 

Traditionally, anisotropy has been quantified as the ra- 
tio (v/a)* of the measured quantity iWx/oo.s to the height 
of the isotropic-rotator curve at a value of e that is char- 
acteristic of the galaxy. This work suggests that anisotropy 
can in future be more rigorously quantified by the anisotropy 
parameter S. 
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